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Abstract. We consider the problem of existence of certain symmetrical solutions of Stokes 
equation on a three-dimensional manifold M with a general metric possessing symmetry. These 
solutions correspond to unidirectional flows. We have been able to determine necessary and 
sufficient conditions for their existence. Symmetric unidirectional flows are fundamental for 
deducing the so-called Darcy's law, which is the law governing fluid flow in a Helc-Shaw cell 
embedded in the environment M. Our main interest is to depart from the usual, flat background 
environment, and consider the possibility of an environment of arbitrary constant curvature K 
in which a cell is embedded. We generalize Darcy's law for particular models of such spaces 
obtained from R 3 with a conformal metric. We employ the calculus of differential forms for a 
simpler and more elegant approach to the problems herein discussed. 
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1. Introduction 

Pattern formation is a very exciting and fastly growing area in physics and related sciences 
[H El El EH Ej The Saffman- Taylor [B] problem is one of the most studied among the systems 
presenting formation and evolution of patterned structures. It studies the hydrodynamic 
instabilities at the interface separating two immiscible fluids confined between two parallel 
flat plates, the Hele-Shaw cell. In such configuration, when a low-viscosity fluid displaces a 
higher viscosity fluid, the interface becomes unstable, deforms, and forms fingers ff\ l8| 19} ITU]. 
The key point in the study of such patterns is Darcy's law, the two-dimensional reduction 
of the Navier-Stokes equation incorporating the boundary conditions (no-slip) and the mass 
conservation law (continuity equation). An interesting generalization of this problem is to 
vary the geometry of the cell and study how geometric parameters like curvature, for instance, 
affect the interface dynamics. Spherical [TT) IT2] , cylindrical ^3] and conical geometries 
have been studied in this context, yielding substantial information connecting relevant fingering 
mechanisms (finger competition and finger tip-splitting) to the cell's geometric and topological 
features. It is our aim here to provide the means of further generalizations of the Saffman- Taylor 
problem involving other geometries. 

In this work we derive Darcy's law in more general geometries, where we allow the ordinary 
three-dimensional space K 3 to acquire a symmetrical Riemannian metric then look for the 
conditions wherein Navier-Stokes equation in this space becomes separable. An important 
case where this occurs is that of a separable metric. For this case, we demonstrate that 
symmetric unidirectional flows are possible. Another important case is that of M 3 endowed with 
a conformal metric. For this case, we derive Darcy's law even in the situations where symmetric 
unidirectional flows do not exist. The case of a pseudo- Riemannian metric (Minkowski) is 
also approached for it supports surfaces of constant negative curvature (pseudospheres or 
Lobachevsky planes), which have been partially studied in [TT]. The calculus of differential 
forms is used throughout the article for a simpler and more elegant way of presenting the 
problem. 

Suppose M is a smooth orientable manifold of dimension 3 endowed with a metric locally 
given by 

ds 2 = E\dx\ + E\dx\ + E\dx\, (1) 

where E^ are smooth functions of the coordinates x^ % — 1, 2, 3. Let ( , ) be the inner product 
induced by this metric. 

The motion of a fluid in M is described by a vector field V : M — > TM, where TM 
denotes the tangent bundle to M. There exists a canonical correspondence between TM 
and the cotangent bundle TM* = /\ 1 (TM*), which is defined using the metric to each 
tangent vector V G TM there corresponds a unique differential 1-form G /\ 1 (TM*) such 
that uj^{W) = (V,W). Given an orthonormal basis f3 = {ei,e 2 ,e 3 } of TM, we define the 



corresponding basis of /\ 1 (TM*) as (3* = {uj Sl) ujg 2) ug. j }. The basis (3* is orthonormal in the 
inner product ( , ) induced by the dual metric 

(ds*) 2 = (l/El)dx\ + {l/El)dx 2 2 + {l/El)dxl. (2) 

Fluid motion is governed by the Navier-Stokes equation 
dV 



P 



—grad(p) + r/AV. 



If the flow is incompressible, V must also satisfy the equation of continuity 
divV = 0. 

In many applications, such as the study of Hele-Shaw flows, one assumes a steady flow and 
neglects the so-called inertial terms on the left-hand-side of the Navier-Stokes equation. Under 
these hypotheses, Navier-Stokes equation reduces to Stokes equation. Using the correspondence 
between vector fields and differential forms, Stokes equation and the equation of continuity 
translate respectively into 

- Ugrad(p) + r}Au$ = ' ( 3 ) 

<W =0, (4) 
where A is the Laplace operator 
-(d5 + 5d). 

We recall that the operators d and 5 are the exterior differential and codifferential, respectively. 
The codifferential is an operator from /\ k (TM*) to /\ k ~ l (TM*) defined by 

5= (-l) k *d*, 

where * : f\ k (TM*) -> f\ 3 ~ k (TM*) is the Hodge star operator. 

— * 

Darcy's law is obtained by averaging V in the normal direction with respect to a given 
two-dimensional smooth submanifold N of M. It provides a reasonable description of the 
fluid motion between two non-intersecting neighboring copies of N. Generally such pair of 
submanifolds is said to form a Hele-Shaw cell. Of particular interest, due to their simplicity, are 
the cells formed by level sets {xi = constant} in a local chart of M. All examples of Hele-Shaw 
cells studied so far (planar, cylindrical, conical, spherical) are formed by such submanifolds. 
One procedure for obtaining Darcy's law consists of considering a one-parameter family of 
velocity fields which corresponds to what we called symmetric unidirectional flows. When such 
family exists and its profile function (denoted in this paper by g) is non-constant, a simple 
method provides a quick deduction of Darcy's law. The existence of symmetric unidirectional 
flows is, therefore, a very important issue. One of our main goals is to attempt to overcome 
the non-existence of such flows in a perturbative way. 

In Section 2 we look for solutions of Stokes equations for a symmetric unidirectional 
flow and analyze under which conditions the equation is separable. In Section 3 we find 
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the conditions for the existence of symmetric flows in curved space and in Section 4 we study 
the solutions of Stokes equation (JHJ) and deduce Darcy's law for the following systems of the 
Hele-Shaw type: 

(1) Two nearby pseudo-spheres in Minkowski's 3-space; 

(2) Two parallel planes in M 3 with a conformal metric. 
Section 5 summarizes our main results and conclusions. 



2. Solution of Stokes equation for a symmetric unidirectional flow 

We study the solutions of Stokes equation in M under the following assumptions 

(Al) The level sets S a = {x± = a} and S a+ b = {xi = a + b} are two smooth (non-intersecting) 
surfaces. The vector field e*i is normal to both S a and S a+ b- 

(A2) The coefficients of the metric of M, do not depend on x 2 . 

(A3) Fluid motion is in the direction of x 3 and the velocity field does not depend on x 2 . We 
refer to such motion as a symmetric unidirectional flow. 

A velocity field satisfying assumption (A3) has the form 

V = V 3 (x u x 3 )e 3 . 

Using the correspondence with 1-forms, we obtain 

u$ = V 3 ug 3 = (V 3 E 3 )dx 3 . 

The equation of continuity then writes 

d(V 3 E 1 E 2 ) 

5u)p = — * d * (V 3 E 3 dx 3 ) = — * d(V 3 Eidxi A E 2 dx 2 ) = — * ( -dx\ A dx 2 A dx 3 ) = 0. 

ox 3 

Thus we have 9<yEl ^. 2 ^ = 0, which implies that 

EiE 2 V 3 = g(x x ) (5) 

for some function g which remains to be determined. 

Using the equation of continuity, the Laplacian of u>^ reduces to — 5du^, i.e. 

fd(V 3 E 3 ) 

— * d * d(V 3 E 3 dx 3 ) = — *d * — dxi A dx 3 

\ OXi 

~ * d ' ~E-E 3 ^r dx '< 



which equals 



d ( E 2 d(V 3 E 3 )\ , , d ( E 2 d(V 3 E 3 , . , 

dxi A dx 2 — — — ( — — ) dx 2 A dx 3 



dxi \E 1 E 3 dxi J dx 3 \E X E 3 dx x 

E x d ( E 2 d(V 3 E 3 )\ dx ^ | E 3 d ( E 2 d(V 3 E 3 ) 



E 2 E 3 dx 3 \EiE 3 dx\ J E 1 E 2 dx\ \EiE 3 dx\ 

Recall that, by definition of the gradient, 

dp dp 

Ugrad(p) =dp= ^—dXi + -— dx 3 . 

OX i ox 3 
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Thus, from the expression for the Laplacian of deduced above, Stokes equation is equivalent 
to the following system of equations 

dp Ex d f E 2 d(V 3 E 3 ) \ 



dx\ E 2 E 3 dx 3 \EiE 3 dx\ 
dp E 3 d ( E 2 d(V 3 E 3 ) 



dx 3 E 1 E 2 dx 1 \E X E 3 dx\ 
The unknown function g[x\) in (jSJ) must be such that the system above is satisfied for some 
smooth function p(x\, x 3 ). If we apply the operator d to both sides of Stokes equation 

—dp + rjAuj-p = 

we obtain, since r] ^ 0, 

dAu^ = 

because d 2 = 0. Hence the (local) existence of p satisfying Stokes equation is guaranteed if g 
solves the equation 

dAuj 9 = -d5d (j^r) = 0, 

or, in coordinates, the equation 

d ( Ex d ( E 2 d {&)\\ d I E 3 d ( E 2 9 (&) 



dx 3 I E 2 E 3 dx 3 \ E\E 3 dx\ J J dx\ I E X E 2 dx\ \ E\E 3 dx\ 

which is a linear ordinary differential equation of third order whose coefficients are functions of 
the metric coefficients. Thus the existence of non-constant solutions of equation (JEJ) depends 
on the metric. From now on we will only consider non-constant solutions of (0) • 

2.1. Separable Stokes Equation 

A great simplification is achieved when the metric is such that Stokes equation is reduced to 
a single separable differential equation. It turns out that for some important examples Stokes 
equation reduces to a separable equation (|7jl. This fact motivates the definition below. 

Definition 2.1 We say that Stokes equation is separable if the system of differential 
equations (Q|) and reduces to a separable equation 

The following definition will also be useful in our discussion. 

Definition 2.2 A function f(x±, x 3 ) is separable if it can be written as a product of a function 
of Xi and a function of x 3 . 



Let us suppose that 

d ( E 2 d ( gE 3 



dx 3 \E 1 E 3 dxi \E ± E 2 
or, equivalently 

E 2 d ( gE 3 



(9) 



E 1 E 3 dx 1 \EiE 2 



C( Xl ). (10) 
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for some function C. 

Multiplying both sides of the equation above by g(x±), and denoting by G, we obtain 

°(f)=r(°fK 

Dividing both sides by G 2 , we get 
which in terms of g is 



EiE 2 J J g E X E 2 g 

The last equation above shows that equation (fTU|) holds as long as the ratio -g^r is the product 
of a function of function of x 3 and a function of x x , x 3 . 

Proposition 2.1 Assume that E\ = E\(x\). Then Stokes equation is separable if and only if 
the ratio J| is separable. 

Proof. If Ei = Ei(xi) and if Stokes equation is separable in the sense of definition 12.11 then 
from equation (jHJ) 

E exp(C)exp (J ^dxi) E 

V ' J H(xi)H(x 3 ) 



EiE 2 g E 2 

for some functions H, H. 

On the other hand, if #■ = H(xi)H(x 3 ) then Q holds and equation flJJ) becomes 



Sari \E1E3 dxi 

After dividing by H(x 3 ) and moving the second term to the right hand side, we have 
I dp vH d f E 2 d (j^\\ =1jH d ( 1 d fgH 



H dx 3 dxi \E1E3dx1 \EiE 2 ) ) dxi \EiH dx\ \ Ei 

Clearly each side of the equation above depends on a single variable. □ 

We now discuss a few examples. 

Example 1: Let M = M 3 with the Euclidean metric in cylindrical coordinates (p,<f),z). We 
have that 

Ei = E 3 = 1, E 2 = p. 

From proposition 12.11 it follows that Stokes equation for a flow in the ^-direction reduces to 
equation (JJJ) 

Op Vdfd/g\\ 



dz pdp \ dp \p 



If, instead, we consider flows in the p-direction (radial flow) which are 0-symmetric, then 
Stokes equation is 

^ = V9". (12) 
pdp 

Finally, for z-symmetric flows in the ^-direction, we have 

* »( I'm V (i3) 



dp \pdp 

Example 2: Let M = R 3 with the Euclidean metric in spherical coordinates (r, 0, 9). In this 

case 

Ei — 1, E 2 — r sin 6 1 , E 3 = r. 
Proposition 12.11 implies that Stokes equation for polar flows reduces to 

S m6^ = rjg". (14) 

Example 3: Our arguments so far apply to semi-Riemannian metrics just as well. Let M = R 3 
with the Minkowski metric 

ds 2 M = —dr 2 + r 2 sinh 2 rdcfi 2 + r 2 dr 2 

in pseudo-spherical coordinates (r, </>, r). The analogous of Proposition 12.11 for a semi- 
Riemannian metric implies that Stokes equation reduces to 
dp 

smhr— = rjg". (15) 

2.2. An example of non-separable Stokes equation: flow in R 3 with a conformal metric 
Suppose now M is R 3 endowed with the metric 
ds 2 = f 2 (x 3 )(dxl + dx\ + dx\). 

It turns out that such choices of M and ds 2 provide examples of 3-spaces of arbitrary 
constant curvature [T3]. Indeed, let K be a non-negative real number. If f(xs) = l/(i^x|+l/4), 
a direct calculations reveals that M has curvature K. For f(xs) = 1/(1 + \^Kx^), we have that 
M = {(xi,X2,xs) G M/xz > 0}, the open half-space, has curvature —K. 

Stokes equations © and (J7J) are 
dp 2g'f 

dx-r 71 —' (16) 

IHf (17) 

If we apply the condition (JBJ for existence of functions g, p satisfying the equations above, 



we obtain 



,/'Y , </' 
' 2 T) 9 = T 
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which is a separable equation. The differential equation above has a non-trivial solution if and 
only if f 3 is constant, i.e, the conformal factor must satisfy a differential equation of the 

form 

By setting y = -jg, we obtain the equation y" = — whose solutions are well-known for any 
values of the constant k\. For each solution y = y(x 3 ) there corresponds an / = y^ 1 ^ and a g 
which is a solution of 

g"' + k ig ' = 0. (18) 

Thus a restriction on the metric needs to be imposed in order to have non-trivial solutions 
of Stokes equation. Unfortunately, this restriction applies to the choices of conformal factors 
which give constant curvature. In particular we have that symmetric unidirectional flows are 
not possible in hyperbolic 3-space. 

3. The existence of symmetric unidirectional flows 

Let us address the fundamental question of existence of solutions of equation (JHJ). 

Recall that the existence of unidirectional flows depends on the existence of a solution 
g = g(xi) of equation (|5ji . As we mentioned before, equation (JHJ) is a third order linear ordinary 
differential equation. It can be put in the form 

g'" - A(xi, x 3 )g" - B{x x , x 3 )g' - C{x u x 3 )g = 0. 

We will show that a differential equation such as this one can only have a solution if the 
coefficients A, B, C do not depend on x 3 . We will need the following 

Lemma 3.1 If the linear differential equation 

g" - A(xi,x 3 )g' - B{x u x 3 )g = 
has a solution, then the coefficients A and B do not depend on x 3 . 
Proof. A solution of the equation in the statement satisfies the linear system 

Ag' + Bg = g" 

A 3 g' + B 3 g = 0, 

where the subscript indicates partial derivative with respect to x 3 . Thus the 2x2 determinant 
A B 
A 3 B 3 

is zero. But from A 3 g' + B 3 g = 0, we have that A 3 is equal to a function /i of X\ times B 3 . 
Hence, we must have A = fiB. So 

g" = B^g' + g). 

Therefore B = B(xi) and A = A(xi), as we wanted to prove. □ 
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Lemma 3.2 If the linear differential equation 

g'" - Afa, x 3 )g" - B(xi, x 3 )g' - C(xi, x 3 )g = 

has a solution, then the coefficients A, B and C do not depend on x 3 . 

Proof. Consider the linear system 

Ag" + Bg f + Cg = g'", (19) 
A 3 g" + B 3 g' + C 3 g = 0. (20) 

From the previous lemma, equation ()20j) implies that there exist functions p, and v of X\ such 
that B 3 = pA 3 and C 3 = vA 3 . Integrating these equations with respect to x 3 , we obtain 

B = fjA + j{xt), C = vA + w(xi). 
Substituting into equation (fTHj) . we obtain 

A(g" + fig' + vg) + ( 7 g' + wg) = g'" . 
Therefore, A is a function of X\ only, and the same must hold for B and C. □ 

Lemma 13.21 provides necessary conditions for the existence of symmetrical unidirectional 
flows on a 3-manifold M with a metric ds 2 whose coefficients are E±, E<i and E 3 . For instance, 
the coefficient of g" in equation (jHJ) is 

Lemma [3.21 says that the right-hand side is a function of x\. Hence, a quick calculation shows 
that -gf§r must be separable. We have proved the 

Proposition 3.1 A necessary condition for the existence of a solution of the differential 
equation (0) is that -gf|r is separable. 

Propositions 12.11 and 13.11 imply the main result of this section. 

Theorem 3.1 Suppose E\ = Ei(x\). If Stokes equation has a solution then ^ is separable. 
Conversely, if ^ is separable, then Stokes equation in the metric with coefficients E x , E 2 and 
E 3 is separable and hence it has a solution. 

Theorem 13.11 implies the impossibility of symmetric unidirectional flows in conical or 
toroidal geometries. 

Example 4 : Let M be an open region of M. 3 endowed with the Euclidean metric in conical 
coordinates (r/, (f>, p) 

ds 2 = drf + [j] cos a + p sin a) 2 d(f) 2 + dp 2 , 

where < a < | is a constant. Since E\ = 1 and 

£3 = 1 

E 2 ri cos a + p sin a 
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is not separable, we have that equation (jSJ) has no solution. 

Example 5: Let M be an open region of 1R 3 endowed with the Euclidean metric in toroidal 
coordinates (r, 6, 9) 



ds 2 = dr 2 + (a + r cos 9) 2 d(j) 2 + r 2 d9 2 , 

where a > is constant. Since E\ = 1 and 

£3 = r 

E 2 a + r cos 9 
is not separable, it follows that (jEJ) has no solution. 

If Ei depends on x 3 , it is possible to have solvable Stokes equations which are not separable, 
as we saw in subsection 2.2. 

Remark 3.2 Theorem AS. 1\ imposes a serious restriction on an argument that has been used to 
deduce Darcy's law for Hele-Shaw systems in curved geometries. In the next section we discuss 
a method which provides a perturbed form of Darcy's law which is valid for the example in 
subsection 2.2, somehow bypassing the obstacle imposed by the non-existence of a symmetric 
unidirectional flow. 

4. Darcy's law 

We now deduce Darcy's law for a separable and a non-separable example of Hele-Shaw systems. 
4-1. Pseudo-Spheres in Minkowski space 

This example illustrates well the method of deduction of Darcy's law for separable systems. 

Recall that x± = r, x 2 = <t> and £3 = r are the pseudo-spherical coordinates of M defined 
in example 3 of subsection 2.1. 

Consider the pseudo-spheres S a = {r = a} and S a +b = {r = a + b} in M. A curve going 
from P G S a to S a +b in the r-direction is given by a path A : [0, 1] — > M, where 

r(t) = r(P)+tb,4>(t) = 0(P),r(t) = r(P). 

We average the function V T (r, r) along the path A. 

In order to solve equation (fTKj) . we set both sides equal to a constant C . If we impose the 

non-slip boundary conditions g(a) = g(a + b) =0, then we must have 

C 

g(r) = - — {r-a){a + b-r). (21) 

Z7] 



The average of V T along A is 

V. 



J x V T ds 



Ix ds ' 

where ds is the element of arc-length of A. We have that 

_ = fjV T (\(t))ibdt = 1 f a+b 9^) dr 
Cibdt frsinhr J a r 
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Using expression (|2Tj) for g, it follows that 

V =-—( (r-a){a + b-r) d \ dp 
2br] \J a r J Or' 

Therefore, Darcy's law for two pseudo-spheres S a , S a+ b in Minkowski's 3-space is 
_ b 2 T { h A 

V T = ^-(gradp) T , (22) 

where T (£) = P(l, 2; 4; -b/a) is a hypergeometric function. For comparison, see the appendix 

of m. 

4-2. Parallel planes in M 3 with a conformal metric 

We return to the example discussed in subsection ([2.2)1 . 

Let us consider the planes S a = {x% = a} and S a +b = {x\ = a + b}. We define a curve 
A : [0, 1] -> M 3 by 

= xi(P)+t6,x 2 (t) = x 2 (P),x 3 (t) =x 3 (P). 

where P G In order to simplify our calculations, we will assume that a = and b > 0. 
The average of a function h : 1R 3 — > M along A is 

/ A hds h(\(t))f(x 3 (P))dt 1 r+» 

We will regard ~~ as an averaging operator, which has the property of being linear with respect 
to functions of x%. Besides, we also have 



If we apply to both sides of equations (JHJ), ffTTTf) and ([T7Jl . we obtain 

F 3 = ^ (24) 

* ' " ^ 7. (26) 



dx 3 \f\ 

If we impose the no-slip boundary conditions g(0) = g(b) = 0, then g' — and hence J^- = 0. 
So, from now on, we will drop equation (|2*3jl . 

As we have seen in subsection 12.21 the function g{x\) must satisfy 

</" + ^' = 0, (27) 

where fci = — |/ 3 ^-pj . Equation (|2*7jl is equivalent to 

g" + hg = k 2 

for some constant k 2 . The solution of this equation satisfying the no-slip boundary conditions 
is 

/ sinha(xi — b) — sinhaxi 



1 + 



V sinh ah 



(28) 
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where k\ = —a 2 . The averages of g and g" are thus 



_ h 
ki 



1 — cosh ab 



ab sinh ab 
Using these formulas, we obtain 



i) \ g" 

Therefore Darcy's law is 

v 3 7 



2r/a 



dp 
dx 3 

sinh ab 



2r]a 



-2k 2 
ab 



sinh ab 



— cosh ab 
sinh ab 



+ 



1 — cosh ab ab 



dp 
dx 3 



+ 



1 — coshftfr ab 



(grad p) 



3- 



If we substitute into the formula above the power series expressions for the hyperbolic sine 
and hyperbolic cosine, we get after some manipulation 



V 3 = J - 

V 



(il-|) + (^-|)^ 2 + --- 



l + 2^ + 2^p + ... 



(gradp) 3 , 



(29) 



where a 2 = |/ 3 ^) • 

The important facts about the formula above are 

(I) It is a generalization of Darcy's law for the flat, Euclidean space. Darcy's law in this 
particular case is obtained by setting / = 1 (and a = 0) in (|29|). 

(II) It is defined for arbitrary conformal factors /, even for the ones for which the corresponding 
symmetric unidirectional flows do not exist! 

In the case of a hyperbolic 3-space of curvature —K, we have 



AK 



(1 + VKx 3 ) 2 

Since x 3 > 0, this expression is less than or equal to 4K, and thus a near-zero choice of 
curvature will make a 2 uniformly small. 



5. Summary and concluding remarks 

In this article we have considered the generalization of fluid flow to non-Euclidean spaces by 
obtaining Stokes equation for symmetric unidirectional flows in a smooth orientable manifold 
of dimension 3. We have also found the conditions under which Stokes equation is separable. 
As examples we recovered Stokes equation in M 3 with the Euclidean metric both in cylindrical 
and spherical coordinates. This was also done for R 3 with Minkowski metric. We studied then 
the case of a flow in M. 3 with a conformal metric and found that a restriction on the conformal 
factor is needed in order to have non-trivial solutions of Stokes equation. This restriction rules 
out spaces of constant curvature such as the hyperbolic 3-space. The existence of symmetric 
unidirectional flows was addressed and a condition on the manifold metric established for Stokes 
equation to have solutions. This conditions rules out symmetric unidirectional flows in conical 
and toroidal geometries. In the conical case, the problem seems to come from the curvature 
singularity at the cone vertex. In reference [Hj this problem was avoided by cutting out the 
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vertex in order to provide the inlet for the flow. Darcy's law was finally obtained for the cases 
of two pseudo-spheres in Minkowiski space and for two parallel planes in IR 3 with conformal 
metric. The latter case considers even the case where symmetric unidirectional flows are not 
possible. A series expansion of Darcy's law for small values of the parameter of separation of 
the two parallel planes recovered Darcy's law in Euclidean space in the unit conformal factor 
limit. Since Darcy's law is the starting point for the study of key nonlinear aspects of the 
Saffman- Taylor problem like finger competition and finger tip-splitting, we hope our work will 
motivate further investigations of such important interfacial features in a variety of curved 
Hele-Shaw geometries. 
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